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O ■ ABSTRACT 

[ We present the first self-consistent theoretical models of magnetar spectra that take 

I into account the combined effects of the stellar atmosphere and its magnetosphere. We 

' find that the proton cyclotron lines that are already weakened by atmospheric effects 

ly-^ ' become indistinguishable from the continuum for moderate scattering optical depths 

■ in the magnetosphere. Furthermore, the hard excess becomes more pronounced due to 
. resonant scattering and the resulting spectra closely resemble the observed magnetar 
, spectra. We argue that while the absence of proton cyclotron lines in the observed 
I spectra are inconclusive about the surface field strengths of magnetars, the continuum 

"fH I carries nearly unique signatures of the field strength and can thus be used to infer this 

Qh| quantity. We fit our theoretical spectra with a phenomenological two-blackbody model 

O ' and compare our findings to source spectra with existing two-blackbody fits. The field 

• strengths that we infer spectroscopically are in remarkable agreement with the values 

■ inferred from the period derivative of the sources assuming a dipole spindown. These 
J> . detailed predictions of line energies and equivalent widths may provide an optimal 

^ , opportunity for measuring directly the surface field strengths of magnetars with future 

^ I X-ray telescopes such as Constellation-X. 



1. Introduction 



Anomalous X-ray Pulsars (AXPs) and Soft Gamma-ray Repeaters ( SGRs) are two similar 



classes of neutron stars thought to be powered by ultras trong magnetic fields ( Duncan Thompson 



2OO2I ). Initially identified by their pulsed soft X-ray (jFahlman &: Gregory I Il98ll ) and 7-ray emis- 
sions (jMazets &: Golenetskii Ill98ll ). respectively, they are by now obse rved in multiple wavebands , 
exhibiting a variety of interesting and complex phenomena (see, e.g., IWoods &: Thompson I l2004l : 
Wang et al. 1120061 : ICamilo et al. 1120061 ). However, their persistent thermal- like X-ray spectra (with 
kT ~ 0.3-0.6 keV) and their X-ray luminosities {Lx ^ 10'^^~^^ e rgs/s) that significantly exc eed their 
spin-down luminosities remain as a defining characteristic (e.g.. iMereghetti &: Stella Ill995l ) and has 
been the subject of numerous studies. 
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Empirical models of the observed X-ray spectra of AXP & SGRs have typically required, in 
addition to a blackbody component, a power-la w component, with photon indices ranging between 
2.0 and 4.5, or a second blackbody (see, e.g., IWoods k. Thompson 1 12004 ) . Although such mod- 
elling is adequate for descriptive purposes, more detailed and physically motivated calculations 
have recently been performed, focusing on the effects of either neutron star atmospheres or their 
magnetospheres on the X-ray spectra. 

Radiative equilibrium models of highly magnetized neutron star atmospheres have been de- 
veloped recently by numerous authors focusing on differen t phys i cal pr ocess es. Fu ll angl e- and 



polarization-dependent scattering has been incorporated bv lOzel I (j200ll ) and iLlovd I (|2003l ). The 
effects of vac uum p o lariza tion resonan c e wer e treated completely, inc luding the r esona nt coupling 
of modes, by lOzel I (120011 ) and iLloyd I (j2003l ). and approximately bv lHo &: Lai I (120031) through a 



probabilistic mo de cor i version app roach. Ion cyclotron lines were calculated by lZane et al. I ()200ll ) 



Ho &: Lai I (|200ll ). and lOzel I (120031). Finally , the effects of the two-dimensional structure of the mag- 



netic field were considered by lLloyd I ([2003) and partial ionization in the atmosphere by lHo &: Lai 



(|2003l ). Most models predict significant deviations from a Planckian spectrum with hard excess that 
depends on the magnetic field strength as well as absorption features at the ion cylotron frequency 
that are somewhat suppressed by the vacuum polarization resonance. 

The surface emission from a neut r on sta r is processed by its magnetosphere before it reaches a 
distant observer. iGoldreich Julian I (119691 ) showed that there is a minimum plasma density given 
by n = 7 X 10~^B/P cm~^ in the magnetosphere of a neutron star (where B is the magnetic field 
strength in Gauss and P is the stellar spin period in secon ds). If the magnetospheric p lasma density 
is equal to this minimum Goldreich & Julian density, iMikhailovskii et al. I ( 19821 ) showed that 



scattering of the surface e mission by such a plasma is not important in the X-ray regime. However, 
in the case of magnetars, [Thompson et al. I (|2002l ) showed that large-scale currents fiowing in the 
magnetosphere result in much larger particle densities such that resonant cyclotron scattering can 
be an effective mechanism to modify the emitted spectrum from the surface of a neutron star. This 
has been proposed as an alternative to the surfa ce models to explain the p ower-law tail of X-ray 
spectra of AXPs and SGRs. Inspired by this idea, iLyutikov &: Gavriil I ([200a) calculated the effects 
of resonant cyclotron scattering on a photon emitted by a Planckian source in the magnetosphere 
of a neutron star. They showed that the initial spectrum can be shifted to higher energies, which 
gives rise to a high-energy tail. 

In the magnetar context, the combined signatures of the atmosphere and the magnetosphere 
on the spectrum have not been studied to date. In this paper we calculate the emitted spectra from 
the surface of a highly magnetized neutron star by incorporating all the physical processes that take 
place in its atmosphere as well as the resonant cyclotron scattering in its magnetosphere. We study 
the effects of scattering, both on the continuum of the spectrum and the proton cyclotron lines. 
We use our results to argue that the absence of the proton cyclotron line is not a good indicator 
of the stellar field strength while the shape of the continuum can be used to infer a combination of 
the effective temperature and the surface magnetic field strength. 
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Models 



We follow the methods discussed in (|6zel I I2OOII . l2003l ) and iLvutikov fc Gavriil I (|2006l ^ to 
calculate the emission from a neutron star atmosphere in radiative equilibrium and the resonant 
cyclotron scattering of photons in the magnetosphere, respectively. 

The neutron star atmosphere models are calculated for fully-ionized H plasmas, taking into 
account the effects of vacuum polarization and ion cy clotron lines. H e re, we also incorporate the 
higher harmonics of the ion cyclotron lines, following iGiiver &: Qzel I (j2006l ). As in the previous 
papers, we construct radiative equilibrium atmospheres using a modified Lucy-Unsold algorithm. 
The surface emission spectrum is completely defined by the effective temperature Tgff of the atmo- 
sphere, the surface magnetic field strength B, and the gravitational acceleration g on the stellar 
surface. Throughout this paper, we set g=1.9xl0^ 



^^'^ cms ^. 



We calculate the effect of the resonant scatt ering on the spectrum using the Green's function 
approach described in iLvutikov &: Gavriil I (j2006l ). We assume that the magnetic field in the mag- 
netosphere is spherically symmetric, following a dependence. We solve the radiative transfer 
equation using a Schwarzschild-Schuster two-stream approximation. The emerging spectrum de- 
pends on two parameters: the resonant scattering optical depth 



2 2 
TT e nr 



(1) 



SmcuJB 

and P, the thermal electron velocity in units of the speed of light. Here, n is the plasma density, 
r is the radius of the cyclotron scattering region for a given photon energy, lob = eB/mc is the 
cyclotron frequency, e and m are the electron charge and mass, respectively. 



3. Results 

In order to investigate the effects of scattering on the surface emission spectrum, we have 
computed model spectra for surface temperatures 0.1-0.5 keV, and for magnetic field strengths 
covering a range from 10^^ to 10^^ G. We have then processed the emission spectrum through the 
magnetosphere using the resonant cyclotron scattering model, for optical depth values r = 1 - 10 
and electron velocity values /? = 0.1 - 0.5. 

To illustrate the effects of resonant cyclotron scattering on both the continuum spectrum and 
the proton cyclotron lines, we plot, in Figure [TJ the spectrum emerging from a neutron star with 
0.3 keV surface temperature, 10^"^ G surface magnetic field strength, and increasing scattering 
optical depth in the magnetosphere. At photon energies above the peak of the spectrum, the non- 
Planckian tails already present in the surface spectra become harder. This is caused by the energy 
gain experienced by the photons at each resonant scattering of the electrons. The deep proton 
cyclotron absorption features present in the surface spectra are smoothed out as a resu lt of the 



effect of resonant cyclotron scattering, as also pointed out by ILvutikov &: Gavriil I (|2006l ). This is 
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Fig. 1. — The spectra emerging from a magnetar with a surface magnetic field strength of 10 G, a 
surface temperature of 0.3 keV, and different scattering optical depths in the magnetosphere. The 
equivalent width of the proton cyclotron line decreases and the spectrum at high energies hardens 
with increasing optical depth. 



- 5 - 



in addition to the fact that vacuum polarization re sonance in th e atmosphere already reduces the 



equivalent widths of the proton cyclotron features iQzel I (|200ll ) . Note that in Figure [H we have 
chosen a surface emission spectrum with particularly strong harmonics of the proton cyclotron line 
to demonstrate the effects discussed here. 

In Figure [2] and Table [H we show quantitatively the reduction of the equivalent width of the 
fundamental proton cyclotron line caused by magnetospheric scattering. The atmosphere model 
calculation in Figure [2] is for a neutron star with a surface magnetic field of 10^^ G and a surface 
temperature of 0.3 keV. Results for different magnetic field strengths and the harmonics of the 
cyclotron lines are presented in Table [TJ As can be seen in Figure [21 the equivalent widths depend 
very weakly on the electron velocity /?, but rapidly decrease with scattering optical depth r. Even 
a scattering optical depth of unity is enough to reduce the equivalent width of the lines almost 
by half. For the model parameters used here, this optical depth corresponds to a particle density 
of ~ 3 X 10^^ cm~^, which is ~ lO"^ times larger than the Goldreich- Julian density for a 10^^ G 
magnet ar with a period of 6s. 

Both vacuum polarization and resonant cyclotron scattering in the magnetosphere have similar 
effects on the X-ray spectrum of a magnetar. These effects give rise to both reduced widths 
of cyclotron lines and to non-planckian shapes in the high energy part of the spectrum (~2.0- 
10.0 keV). Usually this component is modeled phenomenologically by adding a second blackbody 
or a power law to a thermal model. In Table El we show the results of fitting our calculated spectra 
with a phenomenological two-blackbody model. Typically, the temperature of the soft blackbody 
corresponds to the surface temperature of the neutron star, whereas the second blackbody is used 
to fit the high-energy tail of the spectrum. 



Discussion 



The X-ray spectra of magnetars carry signatures of their magnetic field strengths, which deter- 
mine both the shape of their continuum and the energies and equivalent widths of cyclotron lines. 
While the detection of cyclotron lines in AXPs and SGRs would yield the most direct measurement 
of their surface magnetic field strengths, numerous recent searches with Chandra and XMM-Newton 
have typically resulted only in upper limits that are given in Table [3l The only exception has been 



the absorption line feature reported for SGR 0526—66 by llbrahim et al. I (|2003l ). 



Proton cyclotron lines that are more prominent than the current detection limits are predicted 
by all models of the atmospheric emission from magnetars des pite the presence o f processes within 



the stellar atmospheres that reduce their strengths (see, e.g., IZane et al. I l200ll : IHo k. Lai 1 120031 : 



Ozel II2OO3I ). This has been used to argue that the surface field strengths of these sources are such 
that the cyclotron lines are outside of the observed energy ran ge and may be diffe rent than the 



field strengths inferred from the dipole spin-down formula (e.g.. iTiengo et al. II2005I ). 



The theoretical predictions, however, are different when we consider the effects of the mag- 
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Fig. 2. — The change in the equivalent width of the fundamental proton cyclotron line as a function 
of optical depth r, for different values of the electron velocity (3. + sign denotes the equivalent 
width of the line when no scattering taken into account. The atmospheric calculation is for a 
magnetic field strength of 4 x 10^^ Gauss and a surface temperature of 0.3 keV. 
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netosphere. As we showed in §3, the proton cyclotron hnes are smeared and often become indis- 
tinguishable from the continuum even for moderate scattering optical depths. This is because the 
scattering of the surface photons in the resonant layers in the magnetosphere broadens any spectral 
feature due to the random motions of the resonant electrons. It is, therefore, not surprising that 
searches for proton cyclotron lines typically yield no detections. This also shows that the absence 
of cyclotron lines is not a good indicator of the magnetic field strength of the neutron star. 

The magnetic field leaves its imprint also on the continuum via both atmospheric and mag- 
netospheric processes. On the one hand, these two processes work in tandem to strengthen the 
observable continuum signature: Compton upscattering hardens the high energy excess already 
present because of atmospheric effects. On the other hand, the exact shape and hardness of the 
resulting spectrum depends also on the strength of the unscattered cyclotron line, especially if it 
occurs at energies above the peak of the X-ray emission. This is because, even though the line is 
not visible as a distinct feature, it still causes a suppression of the continuum spectrum. As a result, 
when we use a phenomenological two-blackbody spectral model to describe the theoretical spectra, 
the inferred blackbody temperatures depend strongly and not monotonically on the magnetic field 
strength. We show this effect in Figure [3l Note that we have not applied redshift corrections to 
the temperatures obtained by fitting two blackbodies to our theoretical spectra. Applying such 
a correction only shifts the predicted correlations between the temperatures along the diagonal 
without affecting the qualitative result. 

The qualitative dependence of the two blackbody temperatures Ti and T2 changes dramatically 
at magnetic field strength ~ 5 x 10^^ G, where the unscattered cyclotron line is the strongest (see 
Tabled]). For all values of the magnetic field strength, the two temperatures in the phenomeno- 
logical model increase with increasing effective temperature. However, the strength of the second 
blackbody initially decreases with increasing magnetic field, becoming negligible at moderate effec- 
tive temperatures for S ~ 5 x 10^^ G. At even stronger fields, the second blackbody again becomes 
more prominent. 

In Table [21 we show the temperatures of the two blackbodies that represent the observed X-ray 
spectra of several AXPs and SGRs and compare them in Figure 3 to the theoretical predictions. It 
is remarkable that the temperatures of the two blackbodies necessary to fit the source spectra lie in 
the same narrow range predicted by our theoretical models. More importantly, for the two sources 
for which multiple observations requiring different temperature pairs exist, the data nearly follow 
constant field lines. This may provide insight into the evolution of the magnetic field following 
strong bursts. 

Fitting directly the actual theoretical models to the observed spectra of magnetars will provide 
a good handle on their surface magnetic field strengths. However, even the preliminary comparison 
shown in Figure [3] results in magnetic field strengths that are remarkably cl ose to those inferred 



from the spin-down formula for these sources: 2.7 x 10 G for XT E 1810—197 (jGotthelf &: Halpern 



20061 ) ■ (2.4 - 5) X lO^"^ G for IE 1048.1-5947 (jKaspi et al.ll200lh . This implies that very weak (see 
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0.1 0.2 0.3 0.4 0.5 0.6 

(keV) 

Fig. 3. — The temperatures of two black bodies used to fit our theoretical spectra for models with 
different magnetic field strengths in units of 10^^ G (solid lines) as well as the temperatures inferred 
from multiple observations of three magnetars (filled markers). 
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Table 1) cyclotron lines are likely to be present in the X-ray spectra of these sources. Detection of 
such lines with future X-ray telescopes such as Constellation-^ will confirm these predictions with 
direct measurement of the surface magnetic fields of magnetars. 



T. G. wishes to thank to the members of the University of Arizona Physics Department and 
Theoritical Astrophysics Program for their hospitality during the visit to Tucson, Arizona. 
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Table 1. Cyclotron line equivalent widths^ 



Magnetic Field Line Energy Equivalent Width (eV) 



(10^4 G) 


(keV) 


r = 1 


r = 3 


r = 5 


0.1 


0.063 


0.87 


0.44 


.b 


0.1 


0.126 


0.80 


0.41 


.b 


1 


0.635 


27.96 


14.61 


6.7 


1 


1.263 


19.78 


9.54 


3.9 


4 


2.529 


35.05 


12.02 


_b 


4 


5.051 


130.29 


44.07 


14.6 


8 


5.050 


129.68 


27.56 


.b 


10 


6.295 


122.77 


_b 


_b 



^The electron velocity is (3 = 0.3. 

^The lines are indistinguishable from the continuum 
spectra. 
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Table 2. Phenomenological two-blackbody description of magnetar 

models^ 



Magnetic Field 


Surface Temp. 


r 


= 1 


r 


= 3 


T 


= 5 


(10" Q) 














]<T-> 


1 


U.i 


nil 
U.ii 


U.OO 




n /I o 
U.4Z 


U.io 


n Rr\ 
U.oU 


A 


U.i 


U.Uo 


n 1 /I 
U.i4 


U.iU 


n 1 n 

u.iy 


U.iU 


n on 
U.zU 


A 

4 


U.i 


u.uy 


n 1 Q 
U.io 


U.uy 


n oi 
U.Zi 


U.IU 


n OQ 
U.Zo 


/? 
D 


U.i 


U.uy 


n on 
U.zU 


U.iU 


n oo 
U.ZZ 


U.iU 


n oc; 
U.Zo 


Q 

o 


U.i 


n no 

U.uy 


n on 
U.ZU 


U.iU 


n OQ 
U.Zo 


U.iU 


n oc^ 
U.Zo 


1 n 


U.i 


n no 

U.uy 


n on 
U.zU 


n no 

U.uy 


n oo 
U.ZZ 


U.IU 


n O/i 
U.Z4 


i 


U./ 


U.i ^ 


U.ol 


U.io 


A ^; 1 
U.Di 


u.iy 


A ^;a 


o 
Z 


U.Z 


U.ZZ 


U.DD 


U.zo 


A lO 

U./o 


U.ZO 


1 A7 

i.UY 


4 


U.Z 


nil 
U.ii 


n oT 
U.Z ( 


U.i4 


n Qo 
U.oZ 


U.iD 


n Q 
U.OO 


D 


U./ 


n 1 
U.io 


n Qo 
U.oZ 


U.I ( 


n Q7 
U.O / 


U.iy 


n /) 1 
U.4i 


Q 


U.Z 


n 1 7 
U.i / 


n Q 
U.OO 


U.io 


n /in 
U.4U 


U.iy 


n /I /I 
U.44 




U.Z 


n 1 ^7 
U.i / 


n 

U.OO 


U.io 


n /) 1 
U.4i 


U.iy 


U.4o 


1 


U.d 


n oo 
U.ZZ 


n fio 
U.oZ 


U.z4 


n 'TK 
U. / 


U.ZO 


n QQ 
U.OO 


o 
Z 


U.O 


n o 1 

U.ol 


A O O 
U.OO 


n oo 
U.»5z 


A An 

u.yu 


n oo 

U.cJz 


A An 

u.yy 


4 


U.O 


n Q 
U.oO 


_b 


n OQ 
U.OO 


_b 


U.4o 


_b 


/? 
D 


U.d 


U.iD 


n QQ 
U.OO 


U.zi 


U.44 


n o/i 
U.Z4 


n AC\ 

u.4y 


8 


0.3 


0.19 


0.41 


0.23 


0.48 


0.25 


0.53 


10 


0.3 


0.22 


0.44 


0.24 


0.51 


0.26 


0.56 


1 


0.4 


0.25 


0.70 


0.25 


0.80 


0.27 


0.88 


2 


0.4 


0.33 


0.86 


0.35 


0.92 


0.38 


0.99 


4 


0.4 


0.49 


.__b 


0.51 


_b 


0.54 


„b 


6 


0.4 


0.46 


_b 


0.51 


_b 


0.55 


_b 


8 


0.4 


0.46 


_b 


0.27 


0.57 


0.31 


0.61 


10 


0.4 


0.20 


0.48 


0.28 


0.57 


0.31 


0.63 


1 


0.5 


0.26 


0.75 


0.31 


0.87 


0.34 


0.96 


2 


0.5 


0.44 


1.01 


0.44 


1.07 


0.51 


1.23 


4 


0.5 


0.60 


1.32 


0.60 


1.47 


0.63 


1.60 


6 


0.5 


0.53 


b 


0.58 


b 


0.63 


b 


8 


0.5 


0.54 


b 


0.59 


b 


0.65 


b 


10 


0.5 


0.56 


_b 


0.30 


0.63 


0.35 


0.70 



''The electron velocity is /? = 0.3. 
''A second blackbody is not required. 
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Table 3. Observational upper limits on 


line equivalent widths 




Energy Range 


Line Width 




Equivalent width (eV) 




(keV) 


(eV) 


SGR 1806-20^ 


SGR 1900+14^ IE 1048.1-5947^ 


4U 0142+61= 


1 - 2 





- 


1 ^ 


±u 


- 




100 


- 


OU 


on 


15 




200 


- 


OU 




- 


2 - 4 





10 


20 


20 


- 




100 


20 


30 


30 


- 




200 


60 


60 




70 


4 - 5 





20 


20 


30 






100 


30 


35 


55 






200 


40 


50 








400 








200 


5 - 6 





20 


50 


60 






100 


30 


75 


60 






200 


40 


110 








500 








500 


6 - 7 





10 


90 


90 






100 


25 


125 


150 






200 


30 


170 






7- 8 





10 


35 


150 






100 


30 


85 


220 






200 


30 


100 







Mereghetti et al. I ((20061) 



TiengQ et al. I (|2006l l 



Juett et al. I ([2003) 
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Table 4. Inferred temperatures of two-blackbody fits to 
magnet ar spectra 



Source Name 


kTi 


kT2 


Reference 




keV 


keV 




1E1048. 1-5947 


0.47 


1.00 


Tieneo et al. f2005) 




0.44 


0.86 






0.37 


0.76 






0.58 


1.25 


Oosterbroek et al. fl998) 


1E1841-045 


0.47 


1.50 


Morii et al. f2003l 


XTE 1810-197 


0.25 


0.68 


Gotthelf k Habern ^20061 




0.29 


0.71 






0.27 


0.70 






0.25 


0.67 






0.22 


0.60 






0.20 


0.52 






0.19 


0.46 





